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1. Introduction 

Unlike QED, quantum chromodynamics (QCD) forms a color-electric flux-tube between the 
quark and the antiquark in mesonic systems, and this one-dimensional squeezing of the color- 
electric field leads to a linear confinement potential in the infrared region (TJ]. Actually, apart from 
the color-Coulomb energy around quarks, the flux-tube formation has been observed in lattice QCD 
both for QQ [|] and 3Q systems [g, g, |). 

In the flux-tube picture of hadrons, which is idealized as the string picture in the infrared 
region, one can expect "stringy excitations" of hadrons, as shown in Fig.l. This stringy mode 
is non-quark-origin excitation, and therefore it can be regarded as a gluonic excitation. Such a 
gluonic-excited state would be interpreted as hybrid hadrons (qqG and qqqG), which are interesting 
hadrons beyond the quark-model framework [|]]. 

The flux-tube picture of hadrons 
quark ^ ^ anti-quark 

String picture of hadrons 
quark # # anti-quark 



Stringy excitation of hadrons 




Figure 1: Schematic illustration of the stringy excitation of hadrons. The flux-tube picture of hadrons is 
idealized as the string picture in the infrared region, which is expected to allow "stringy excitations" of 
hadrons. Since this stringy mode is non-quark-origin excitation, it can be regarded as a gluonic excitation. 



In lattice QCD, from a detailed Wilson-loop analysis, the excited-state potentials and the glu- 
onic excitation have been calculated both for spatially-fixed QQ systems and for 3 Q systems 
[j7|, §]. For simpler QQ systems, the behavior of the gluonic excitation is almost consistent with the 
string excitation in infrared region, in spite of a significant difference at the small distance 

In the previous work, IR/UV-gluon contribution to the ground-state potential has been studied 
[Q, and the confinement force is found to be almost unchanged even after the cut of high- 
momentum gluon components above 1.5GeV [^, [K]] in the Landau gauge. This means that the 
confinement property is insensitive to UV gluons. 

In this work, we study not only ground-state potential but also low-lying even-parity excited- 
state potentials of QQ systems in terms of gluon momentum component in the Coulomb gauge. By 
introducing UV-cut in three-dimensional gluon-momentum space, we study the UV-gluon contri- 
bution to excited-state potentials and stringy excitations. 
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2. Lattice formulation 

2.1 Formalism to extract excited-state potentials in lattice QCD 

We present the formalism to extract the excited- state potential ^ ||] for the spatially-fixed QQ 
system. We denote the nth eigen- state of the QCD Hamiltonian Hby \n), 

H\n) = V n \n). (2.1) 

Here, V n denotes nth excited-state potential, and Oth eigen-state means the ground-state. Consider 
arbitrary independent QQ states \(j>k) (k = 0, 1,2...). Generally, each QQ state can be expressed 
by a linear combination of the QQ physical eigen-states: 

|^)=4|0)+^|1)+4|2) + ... . (2.2) 

The Euclidean-time evolution of the QQ state \<j>k(t)) is expressed with the operator e~ Ht , 
which corresponds to the transfer matrix in lattice QCD. The overlap (0 7 (r)|0^(O)) is given by the 
Wilson loop Wj< k , sandwiched by initial state <pk at / = and final state <pj at / = T, and is expressed 
in the Euclidean Heisenberg picture as 

W* = {<l>j(T)\(l> k (0)) = {<l>j\W(T)\(j> k ) = ((j>j\e- HT \(j>^ (2.3) 

OO CO CO 

= E Zc J m c k n (m\e- HT \n) = £ cj^c*, (2.4) 

m—0n—0 n—0 

with the complex-conjugate notation of c J n = (c J n )*. This is a basic relation between Wilson loops 
and potentials. By introducing the matrices C and A T such that 

c nk = c k^ A mn = e ~V n T § mn ? (25) 

this relation can be rewritten as 

W T = C f A T C. (2.6) 

In general, C is not a unitary matrix, and depends on the choice of \<j>k). Using this relation, 
we extract the potentials V n (n = 0, 1,2- • • ) from the Wilson loop Wt. Consider the following 
combination: 

Wj l W T+ \ = {C^A T Cy l C^A T+l C = C- 1 diag(e- Vo ,e- v \e- V2 ,...)C. (2.7) 

Then, e~ Vn can be obtained as the eigen- values of the matrix Wj 1 Wt+\. In fact, they are the 
solutions of the secular equation, 

detjWf l W T+ i-t\) = Y\( e - Vn -t) =0. (2.8) 

n 

In this way, the potentials V n (n = 0, 1 , 2, ...) can be obtained from the Wilson loop matrix, Wf 1 Wt+\ • 
In the practical calculation, we prepare gauge-invariant QQ states |<^) composed by fat-links 
obtained with APE smearing method [pT|], and calculate many Wilson loops sandwiched by various 
combination of initial state \<f>k) and final state . By solving the secular equation Eq. ( |2lf ) 
within a truncated dimension, ground-state and excited-state potentials can be obtained. 
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2.2 Discrete Fourier transformation and UV-cut of gluon momentum components 

In this subsection, we consider the three-dimensional Fourier transformation of the link- variable 
Un(x) E SU(3) on a periodic lattice of size L 4 , and introduce UV-cut in three-dimensional momen- 
tum space [j^ [lOj] . For the argument on the gluon momentum, gauge fixing is generally needed. 
For the comparison with continuum QCD, the suitable gauge to be taken on lattice would be the 
Landau or the Coulomb gauge, where the gauge field tends to be continuous. 

Here, we consider link- variables fixed in the Coulomb gauge, because spatial gauge-field fluc- 
tuation is strongly suppressed. The Coulomb gauge has a global definition to minimize the "total 
amount of the spatial gauge-field fluctuation", 

R = Jd 3 xtr{Ai(x,t)Ai(x,t)} = i Jd 3 x Af(x,t)Af(x,t) (2.9) 

The Coulomb gauge has a physical meaning that it maximally suppresses artificial fluctuation stem- 
ming from gauge degrees of freedom for spatial gluons. In lattice QCD, the Coulomb gauge fixing 
is expressed in terms of link- variable and is defined by the maximization of 

fl latt = ££Retr^,0. (2.10) 

x i 




Figure 2: A schematic figure of UV-cut with Auv and IR-cut with Air on momentum- space lattice, of which 
spacing is given by a p = 2n/ (La). The momentum- space link- variable U^(p) is replaced by the free-field 
link- variable U f * eQ (p) = S p o in the shaded cut regions. 



Now, we perform the three-dimensional discrete Fourier transformation of the link-variable 
Ufi(x) E SU(3), and define the "momentum- space link- variable": 

Uli(p>t) = l^(*>') ex P ( 2 - n ) 

X 

We introduce "UV-cut" in the momentum space. Outside the cut, (/?, t) is replaced by 0, since 
momentum- space link- variable is U f * Qe (p) = S p o in the free-field case of U f * eQ (x) = 1. (See Fig.l.) 
We define "UV-cut momentum- space link- variable": 

~a/^ x \U u (p,t) for \p\ < Auv 

U£(?,t) = { J (2.12) 

10 for \p\ > A uv 
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By the three-dimensional inverse Fourier transformation 

K (x, t) = £ G£(p, t) exp (-ip. x) 



(2.13) 



and SU(3) projection by maximizing 

Re tr{t#(2, t)U*{%t)}, 
we obtain "UV-cut (coordinate- space) link- variable": 

f) GSU(3). 

Using the UV-cut link- variable Ufr(x) instead of U tI (x), we calculate many Wilson loops Wj sand- 
wiched by various combination of initial state \(j>k) and final state Note here that the UV-cut 
should be introduced also to t/ 4 (jt). Otherwise, the QCD Hamiltonian is not changed, so that the 
potentials V n are not changed at all. 



(2.14) 



(2.15) 



3. Ground-state and excited-state QQ potentials and gluonic excitation energy 

In this section, we show the lattice QCD results of the ground- state/excited- state potentials 
and gluonic excitation energy in QQ systems with or without UV-cut. Our numerical simulation is 
performed with isotropic plaquette gauge action with j8 = 6.0 at the quenched level. The lattice size 
is 16 4 , and the periodic boundary condition is imposed. This lattice QCD condition corresponds 
to the (coordinate- space) lattice spacing a ~ 0.104fm, and the momentum- space lattice spacing 
a p = 2n/ {La) ~ 0.74GeV. We use 100 gauge configurations, and average all the parallel-translated 
Wilson loops in each configuration. 
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Figure 3: Ground-sate and even-parity excited-state potentials in QQ systems with or without UV-cut plotted 
against the interquark distance R. The left panel shows the result without UV-cut. The middle and right 
panels show the results with the UV-cut of Auv = 3<2p ~ 2.2GeV and Auv = 2<2 P c± 1.5GeV, respectively. 
The circle symbol denotes the ground- state potential. The square and the triangle symbols denote the even- 
parity excited- state potentials. 

For simplicity, we only consider even-parity exited-state potentials in this paper. We prepare 
the QQ state 1 0*) (k = 0, 1 , 2, 3) composed by the "fat-links" obtained by the APE smearing method 
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QlTp with the smearing parameter a = 2.3 and the iteration number of N smi = 0, 8, 16, 24. Note that 
only even-parity components can be obtained by this parity-invariant procedure. 



3.1 Ground-state and excited-state QQ potentials with UV cut 

Figure [3] shows ground-state and excited-state potentials in QQ systems with or without UV- 
cut of gluon fields. In the original no UV-cut case, the IR slopes of ground-state and excited-state 
potentials are almost the same, as was indicated by the previous lattice studies This means the 
same confinement force in the infrared region. 

By the cut of UV-gluon above Auv = 3a p ~ 2.2GeV, the short-distance Coulomb part pro- 
portional to 1/r reduces in ground-state potential. In the case of Auv = 2a p ~ 1.5GeV, the short- 
distance Coulomb part disappears in the ground-state potential. These tendencies are consistent 
with the previous studies ^ [l(J. 

On the other hand, the shape of the excited-state potential is largely changed by the UV-cut of 
gluon fields for Auv = 1-5, 2.2GeV, while the ground-state potential is not so changed except for 
the short distance. As a caution, the physical size of our lattice is rather small, and the true IR slope 
of the excited-state is expected to be unchanged, because no change is found in the ground-state 
potential, which gives a lower bound of V n in the infrared region. In any case, the change of the 
excited-potential is more significant than that of the ground-state potential by the UV-cut of gluons. 

In the string picture, this result seems to be natural as mentioned below. For the stringy-excited 
state as shown in Fig.l, there is a typical wavelength proportional to the interquark distance R, and 
this wavelength is smaller for higher excitation mode. Then, we expect a significant influence of 
the removal of UV-gluons for the stringy mode, when the UV-cut length 1/Auv becomes larger 
than the typical wavelength of the stringy excitation. In fact, the effect of UV-gluon cut would be 
larger for higher excitation. Our lattice QCD results seem to be qualitatively consistent with this 
tendency. 



3.2 Gluonic excitation energy with UV cut 

The gluonic excitation energy defined by the difference, V n — Vo, is shown in Fig|| Roughly, 
even after the removal of UV-gluons, the magnitude of gluonic excitation is approximately un- 
changed, and the low-lying gluonic excitation remains to be of the order of lGeV 
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Gluonic Excitation : max03-min00 
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Figure 4: The gluonic excitation energy defined by the difference between the ground- state and the excited- 
state QQ potentials, AE n = V n — Vo. The left panel shows no UV-cut case. The middle and right panels show 
the results with the UV-cut of Auv = 3a p ~ 2.2GeV and Auv = 2<2 P ~ 1.5GeV, respectively. 
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4. Summary and concluding remarks 

In this study, we have studied ground-state and low-lying even-parity excited-state potentials 
of quark-antiquark systems in terms of the gluon momentum component in the Coulomb gauge 
using SU(3) quenched lattice QCD. By introducing UV-cut in the gluon momentum space, we 
study the "UV-gluon sensitivity" of the potentials and the stringy excitation. Even after cutting 
off high-momentum gluon component above 1.5GeV, the IR part of the ground-state potential is 
almost unchanged. On the other hand, the change of excited-state potential is more significant by 
the cut of UV-gluons. However, the magnitude of the low-lying gluonic excitation remains to be of 
the order of lGeV after the removal of UV-gluons. 

As a next step, we will investigate the odd-parity excited-state potential, by using non-symmetric 
state operators of\(j)) k . In this work, we use the Coulomb gauge fixing, however, to remove gauge 
artifact completely, it would be also interesting to apply the gauge-invariant expansion in terms of 
the Dirac mode [0. 
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